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ABSTRACT 

Th  i  s  report  d  escr  i  bes  the  ap  p  I  i  cati  on  of  the  F I  u  ent  cod  e  to  the  nu  meri  cal  si  mu  I  ati  on  of  the  f  ree- 
surfaceflow  around  a  model  naval  ship;  the  DTM  B  5415.  Simulations  were  performed  using  both 
a  structured  hexahedral  mesh  and  an  unstructured  tetrahedral  mesh  of  lower  resolution.  The 
results  show  that  Fluent  is  ableto  accurately  si  mu  I  ate  the  total  ship  resistance,  near-field  wave 
shapes,  and  the  velocity  field  in  the  propeller  plane.  The  results  indicate  that  Fluent  is  a  viable 
tool  which  could  be  considered  for  use  in  more  demanding  naval  problems,  such  as  the 
computation  of  ship  wakes  undergoing  specified  manoeuvres.  A  number  of  significant  problems 
were  encountered  duri  ng  the  course  of  the  work  and  have  been  summarized  to  faci  I  itate future 
applications  of  the  code  to  similar  naval  problems. 
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Fluent  Code  Simulation  of  Flow  around  a  Naval 
Hull:  the  DTMB5415 


Executi  ve  Su  mmary 

The  DSTO  has  a  requirement  to  calculate  propeller  loadings,  performance  and 
signatures  for  both  naval  surface  ships  and  submarines.  A  detailed  knowledge  of 
the  fluid  velocity  in  the  propeller  plane  is  required  in  order  to  perform  these 
calculations.  Modern  Computational  Fluid  Dynamics  (CFD)  computer  codes  which 
solve  the  Reynolds-Averaged  Navier-Stokes  (RANS)  equations  for  complex 
geometries  have  been  used  to  simulate  viscous  flow  around  ship  hulls  si  nee  the 
early  1990s. 

This  report  d  escr  i  bes  the  appl  i  cati  on  of  the  F I  uent  code  to  the  nu  meri  cal  si  mu  I  ati  on 
of  the free-surface flow  around  a  model  naval  ship;  the  DTM  B  5415.  Thefocus  of 
the  work  is  to  study  the  capability  of  the  Fluent  code  to  accurately  simulate  the 
total  ship  resistance,  near-field  waveshapes,  and  the  velocity  field  in  the  propeller 
plane. 

Simulations  were  performed  using  both  a  structured  hexahedral  mesh  and  an 
unstructured  tetrahedral  mesh  of  lower  resolution.  Our  results  show  that  Fluent  is 
ableto  accurately  si  mu  I  ate  the  total  ship  resistance,  near-field  waveshapes,  and 
the  vel  oci  ty  f  i  el  d  i  n  the  propel  I  er  pi  ane.  The  total  resi  stance  coeff  i  ci  ent  cal  cu  I  ated 
on  the  hexahedral  mesh  wasfound  to  agree  with  theexperi  mental  valueto  within 
3. 8%.  The  si  mu  I  ated  waveshapes  along  thesurfaceofthehull  and  in  the  near  field 
computed  on  both  the  hexahedral  and  tetrahedral  meshes  were  found  to  be  in 
good  qual  itative agreement  with  theexperi  mental  profi  les,  with  thedegreeof  error 
from  the  resu  I  ts  on  the  tetrahed  ral  mesh  bei  ng  si  mi  I  ar  to  those  from  the  hexahed  ral 
mesh.  The  axial  velocity  component  in  the  propeller  plane  was  found  to  be  the 
most  difficult  aspect  of  the  flow  to  calculate  accurately,  although  the  maximum 
error  was  never  more  than  10%. 

Theresults  indicate  thattheFluentcodeisaviabletool  which  could  be  considered 
for  use  in  moredemanding  naval  problems,  such  as  the  computation  of  ship  wakes 
undergoing  specified  manoeuvres.  A  number  of  significant  problems  were 
encou  ntered  d u ri  ng  the  cou  rse  of  thi  s  w or k  and  they  are  d escr i  bed  i  n  the  report  to 
facilitate  future  applications  of  thiscodeto  similar  naval  problems. 
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1.  Introduction 

TheDSTO  has  a  requirement  to  calculate  propeller  loadings,  performance  and  signaturesfor 
both  naval  surface  ships  and  submarines.  A  detailed  knowledge  of  the  fluid  velocity  in  the 
propeller  plane  is  required  in  order  to  perform  these  calculations.  Modern  Computational 
Fluid  Dynamics  (CFD)  computer  codes  which  solve  the  Reynolds-Averaged  N  avi  er- Stokes 
(RANS)  equations  for  complex  geometries  have  been  used  to  si  mu  I  ate  viscous  flow  around 
shi  p  hul  I  s  si  ncetheearl  y  1990s.  Thefi  rst  si  gnifi  cant  si  mul  ati  ons  wereconfi  ned  to  tanker  hulls 
and  were  relatively  simple  by  today's  standards;  wave  effects  were  absent,  the  free  surface 
was  considered  as  a  plane  of  symmetry,  and  the  most  sophisticated  turbulence  model 
employed  was  the  k-s  model.  A  summary  of  these  early  simulations  can  be  found  in  the 
proceedings  of  the  SSPA-CTFI -I  I  FIR  workshop  on  Ship  Viscous  Flow  held  in  Gothenburg, 
Sweden,  in  1990 [1]. 

Asthe  speed  of  computers  increased  and  more  sophisticated  RANS  codes  were  developed 
more  real  i  sti  c  si  mu  I  ati  ons  w  ere  abl  e to  be  performed .  These  advances  are  w  el  I  d  ocu  mented  i  n 
the  proceedings  of  several  international  conferences  on  the  application  of  CFD  techniques  to 
ship  flows  which  have  been  held  every  few  years  since  1990,  most  notably  those  in  Tokyo  in 
1994  [2],  Gothenburg  in  2000  [3]  and  Tokyo  in  2005  [4],  Significant  information  on  the 
application  of  CFD  codes  to  naval  ships  and  submarines  can  also  be  found  in  recent 
proceedings  of  the  "Symposium  on  Naval  Flydrodynamics",  which  isabiennial  symposium 
sponsored  by  the  US  Office  of  N  aval  Research  (ON  R)  and  was  first  held  in  1956. 

An  example  of  the  sophistication  of  current  specialised  CFD  codes  for  the  si  mul  ati  on  of  flow 
around  naval  hulls  is  provided  by  the  work  of  Burg  etal.  [5]  reported  atthe24th  Symposium 
on  Naval  Flydrodynamics.  They  used  U2NCLE,  a  three-dimensional  unstructured, 
parallelized  CFD  code  developed  by  the  Computational  Simulation  and  Design  Centre  at 
Mississippi  StateUniversitytosimulatethefreesurfaceflow  around  a  fully-appended  model 
of  a  U  S  naval  shi  p.  The  si  mu  I  ati  ons  prod  uced  by  thi  s  code  real  i  sti  cal  I  y  captu  red  the  turbul  ent 
flow  and  vortices  arising  from  the  bulbous  bow  and  the  tips  of  the  propul  sors  and  rudders 
and  were  found  to  accurately  model  the  rotation  of  the  propul  sors.  Near-field  waveshapes 
were  accurately  simulated  using  a  nonlinear  free  surf  ace  algorithm. 

CDFShip-lowa,  which  was  developed  at  The  University  of  Iowa's  Institute  of  FH ydraul ic 
Research  (I  I  FIR)  underfunding  provided  by  ONR,  is  another  example  of  aRANScodewhich 
has  been  specifi  cal  I  y  devel  oped  for  surface-shi  p  and  mari  ne-propul  sor  fl  ow  probl  ems.  Typi  cal 
applications  in  naval  hydrodynamics  include  the  prediction  of  resistance  (friction  and 
pressure  drag),  wave  profiles,  sea  keeping  (the  code  is  six-degree-of-freedom  capable),  and 
manoeuvring. 

Wilson  et  al.  [6]  recently  compared  the  performance  of  CFDShip-lowa  with  that  of  two 
commercial  CDF  codes,  Fluent,  developed  by  Fluent  Inc.,  and  Comet,  developed  by 
CD-Adapco,  in  the  prediction  of  ship  generated  wave  fields.  It  was  found  that  each  of  the 
codes  had  different  advantages  and  disadvantages,  and  that  each  had  certain  specific 
requi  rements  for  obtai  ni  ng  accurate  sol  uti  ons  of  a  surface  shi  p  wave  fi  el  d .  1 1  was  al  so  noted 
that  the  commercial  solvers  may  have  an  advantage  when  complex  hull  shapes  involving 
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appendages  and  propul sors  are  considered  because  they  have  the  flexibility  to  use 
unstructured  aswell  as  hybrid  meshes.  Both  Fluent  and  Comet  also  have  the  advantage  of 
al  I  ow i  ng  sol  uti  on  based  gri  d  adapti  on  techni  ques to  provi  defi  ner  gri  d  resol  uti  on  i  n  restri  cted 
regions  of  the  grid,  such  as  in  the  air-water  interface  region,  and  hence  may  offer  a  more 
robust  and  computationally  economical  way  to  provide  accurate  free  surface  predictions  in 
the  vicinity  of  surface  ship  hulls. 

The  Hydrodynamics  Group  within  M  PD  has  considerable  experience  in  the  use  of  the  Fluent 
code  and  its  application  to  a  variety  of  flow  related  problems  in  the  naval  context.  Recent 
examples  involving  submarine  related  flows  indudetheevaluati  on  of  hull  modificationsfor 
the  Col  I  i  ns  C I  ass  Submari  ne  [  7,8],  the  anal  ysi  s  of  fl  ow  d  i  sturbances  created  by  fi  n-mounted 
cameras  on  a  Col  I  i  ns  C I  ass  su  bmari  ne  [9]  and  the  si  mu  I  ati  on  of  the  hyd  rodynami  c  forces  on  a 
model  submari  ne  [10],  Moregeneral  applicationshaveindudedthesimulationofflow  around 
a  generi  c  f i  n- body  j  u  ncti  on  [  11]  as  w el  I  as  tu rbu I  ence  i  n  the  w ake  of  a  generi  c  bl  uff  body  f I  ow 
[12]. 

This  report  describes  the  application  of  the  Fluent  code  to  the  numerical  simulation  of  the 
free-surf  ace  flow  around  a  model  naval  ship;  the  DTM  B  5415.  This  model  was  conceived  by 
the  US  David  Taylor  Model  Basin  (DTMB)  in  the  early  1980s  as  a  preliminary  design  for  a 
surface  combatant  with  a  sonar  dome  bow  and  transom  stern.  The  model  has  been  studied 
extensi  vel  y ,  both  experi  mental  I  y  and  nu  meri  cal  I  y,  and  i  s  one  of  the  bench  mark  mod  el  s  for  the 
ship  hydrodynamics  community,  having  been  used  for  software  validation  at  the  Ship 
Hydrodynamics  CFD  Workshops  in  Gothenburg  in  2000  [3]  and  Tokyo  in  2005  [4],  Our  focus 
in  this  work  is  to  study  the  capability  of  the  Fluent  code  to  accurately  simulate  the  total  ship 
resistance,  near-field  waveshapes,  and  the  velocity  field  in  the  propeller  plane. 


2.  Computational  Approaches  to  Free  Surface  Modelling 

An  important  aspect  of  any  simulation  of  the  fluid  flow  around  asurfaceship  isthe  method 
used  to  model  the  air/  water  interface.  There  are  basically  two  approaches  to  this  problem; 
surface  fitting  approaches  or  surface  capturing  approaches.  In  the  surface  fitting  approach 
only  the  water  side  of  the  domain  is  simulated  and  the  grid  is  adjusted  to  conform  to  the 
position  of  the  free  surface.  In  the  surface  capturing  approach,  the  computational  domain 
includes  both  the  water  and  air,  and  the  location  of  the  i  nterface  on  the  mesh  is  computed 
from  the  sol  uti  on  of  an  auxiliary  equation. 

Most  of  the  simulations  in  ship  hydrodynamics  prior  to  2000  used  surf  ace  fitting  approaches 
[3].  Exampl es  of  thi s  type  of  method  i  nd  ude  the  si  mul ati  ons  of  Burg  and  M  arcum  [13],  who 
used  a  nonlinear  free  surface  algorithm  in  an  unstructured  finite  volume  RANS  code  to 
si  mul  ate  flow  around  the  DTMB  5415,  aswell  as  the  work  of  Li  [14],  who  used  a  nonlinear 
free  surface  algorithm  on  a  structured  mesh  to  model  flow  around  a  container  ship,  a  tanker, 
and  al  so  the  DTM  B  5415.  The  advantage  of  these  methods  is  that  the  ai  r/  water  i  nterface  can 
be  defined  with  great  precision.  There  are  significant  disadvantages  however.  In  time- 
dependent  si  mul  ati  ons  with  breaki  ng  waves  the  deformati  on  of  the  i  nterface  is  so  great  that 
changes  in  grid  topology  occur  and  these  are  difficult  to  handle  with  a  surface  fitting 
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approach.  Even  without  a  change  in  grid  topology  the  deformation  of  the  i  nterf  ace  can  be  so 
great  that  strong  cell  deformation  can  occur  and  cause  the  computation  to  fail.  Because  of 
these  disadvantages  there  has  recently  been  an  increase  in  the  use  of  surface  capturing 
techniques  for  the  simulation  of  steady  and  unsteady  ship  motion  in  the  field  of  naval 
hydrodynamics. 

Surface  capturing  approaches  can  declassified  aseitherVolume-of-Fluid  methods(VOF)  [15] 
or  Level  Set  (LS)  methods  [16,17].  TheVOF  method  can  si  mu  I  ate  the  motion  of  two  or  more 
immiscible  fluids  by  solving  a  single  set  of  momentum  equations  and  tracking  the  volume 
fractions  of  each  of  thefluids  throughout  the  domain.  A  volume  fraction  function  is  defined 
for  each  of  thefluids  in  the  simulation  and  is  set  to  unity  if  the  given  fluid  occupies  a  cell 
volume  or  is  set  to  zero  otherwise.  When  the  interface  between  two  fluids  cuts  through  a 
computational  cell  thevalueofthefunction  for  a  particular  fluid  represents  thefracti  on  of  the 
cell  volume  occupied  by  that  fluid.  As  the  simulation  evolves  the  volume  fraction  functions 
are  con vected  by  the  u  nd  erl  y i  ng  f  I  u  i  d  f I  ow .  The  V 0  F  method  therefore  has  the  abi  I  i  ty  to  treat 
quite  complex  interface  evolutions  and  also  has  good  mass  conservation  properties.  The 
method  is  not  ideally  suited  to  problems  involving  surface  tension  effects  however  as 
problems  can  arise  in  constructing  sufficiently  accurate  expressions  for  the  surface  curvature 
and  su rface  normal  vector  from  the  vol  umefracti  on  fundi  ons.  Thi  s  i  s  not  normal  I y  a  probl  em 
in  surface  ship  simulations  however.  Recent  examples  of  the  use  of  the  VOF  method  to 
successf  u  1 1  y  si  mu  I  ate  the  f  I  ow  arou  nd  the  DTM  B  5415  i  nd  u  d  e  the  w  ork  of  C  hen  et  al .  [  18]  and 
Rheeand  Skinner  [19]. 

I  n  the  LS  method  the  i  nterf  ace  between  thetwofluidsis  represented  by  the  set  of  al  I  poi  nts  of 
some  higher  order  fundion  for  which  the  fundion  has  the  value  zero.  This  embedding 
fundion  is  then  evolved  with  the  flow  and  the  zero  level  set  at  any  instant  signifies  the 
location  of  thesurface.  Sincea  smooth  fundion  can  be  chosen  to  represent  the  initial  surface 
I  ocati  on  the  cal  cu  I  ati  on  of  the  cu  rvatu  re  and  other  geometri  c  featu  res  of  the  su  rface  as  thef  I  ow 
evolves  is  relatively  straightforward.  One  disadvantage  of  the  method  however  is  that  it  not 
good  at  conserving  the  mass  of  individual  fluids  on  the  mesh  and  this  can  lead  to  problems 
for  lengthy  time-dependent  simulations. 

I  n  a  two-phase  LS  method  both  ai  r  and  water  aresi  mul  ated  as  a  si  ngl  efl  ui  d  whose  properti  es 
vary  conti  nuously  across  thefree  surface,  whereas  i  n  a  si  ngl  e-phase  LS  method  thesol  uti  on  i  s 
computed  only  in  the  liquid  phase.  There  has  been  significantdevelopmentintheapplication 
of  LS  methods  to  naval  hydrodynamic  problems  in  the  last  few  years.  Yang  et  al.  [20]  have 
implemented  a  two-phase  level  set  method  with  an  immersed-boundary  Cartesian  grid 
method  in  version  six  of  the  CFDShip-lowa  code  and  used  it  to  perform  large-eddy 
si  mul  ati  ons  (L  ES)  of  several  shi  p  geometri  es,  i  nd  ud  i  ng  the  DTM  B-5415,  and  promi  si  ng  results 
have  been  obtai  ned .  D  i  M  asci  o  et  al .  [21]  used  a  si  ngl  e-phase  LS  method  to  si  mul  ate  both  non- 
breaking  and  breaking  flows  around  the  DTM  B  5415  and  Carrica  et  al.  [22]  used  a  single- 
phaseLSmethod  in  combination  with  asixdegree  of  freedom  (6DOF)  algorithm  and  dynamic 
overset  grids  to  calculate  the  sinkage  and  trim  and  pitch  and  heave  motions  for  the  DTM  B 
5512,  which  is  a  smaller  scale  model  of  theDTMB  5415. 

There  are  two  approaches  availablefor  the  numerical  simulation  of  multi  phase  flows  in  the 
Fluent  code;  the Euler-Lagrangeapproach  and  the Euler-Euler  approach.  The Euler-Lagrange 
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ap  p  roach  i  s  essenti  allyaLagrangiandi  screte  phase  mod  el  i  n  w  h  i  ch  the  f  I  u  i  d  p  h  ase  i  s  treated 
as  a  continuum  by  solving  the  N  avier- Stokes  equations,  while  the  dispersed  phase  is  solved 
by  tracking  a  I  argenumber  of  particles,  bubbles,  or  droplets  through  thecal  culated  flow  field. 
This  is  obviously  relevantto  particle-laden  flows  and  isinappropriateforthesea/  air  interface 
problem  considered  here,  although  it  may  become  relevant  in  the  future  if  breaking  waves 
need  to  be  considered.  An  example  of  this  type  of  application  would  occur  during  the 
simulation  of  bubbly  ship  wakes  for  the  calculation  of  torpedo  tracking  signatures. 

In  the  Euler-Euler  approach  the  different  phases  are  treated  as  conti  nua  and  three  different 
Euler-Euler  multiphase  models  are  available;  the  mixture  model,  the  Eulerian  model,  the 
volumeof  fluid  (VOF)  model.  Themixture  model  is  designed  for  two  or  more  phases  in  which 
one  of  the  phases  represents  a  parti  cul  ate.  1 1  solves  for  the  mixture  momentum  equati  on  and 
then  prescribes  relative  velocities  to  describe  the  dispersed  phases.  Applications  of  the 
mixturemodel  include  particle-laden  flows,  sedimentation,  and  bubbly  flows.  The  Eulerian 
model  i  s  the  most  compl  ex  of  the  multi  phase  model  s  i  n  that  it  sol  ves  a  set  of  momentum  and 
continuity  equations  for  each  phase  and  coupling  isachieved  through  pressureand  interphase 
exchange  coefficients  which  depend  upon  thetypeof  phases  involved.  This  model  issuitable 
for  the  simulation  of  bubble  columns,  particle  suspensions  and  fluidized  beds. 

TheVOF  model  isasurface-trackingtechniqueappliedtoafixed  Eulerian  mesh.  Itisdesigned 
for  tw o  or  more  i  mmi  sci  bl  e  f I  u i  d s  w  here  the  posi  ti  on  of  the  i  nterf  ace  betw een  the  f I  u  i  d s  i  s  of 
interest.  A  singleset  of  momentum  equations  is  shared  by  thefluidsand  the  volume  fraction 
of  each  of  the  f  I  u  i  d  s  i  n  each  computati  onal  cel  I  i  s  tracked  th  rou  ghout  the  d  omai  n .  Th  i  s  mod  el 
is  ideally  suited  to  applications  involving  free-surf  ace  flows,  filling,  sloshing  ,  the  motion  of 
large  bubbles  in  a  liquid,  and  the  motion  of  liquid  after  a  dam  break. 


3.  Previous  CFD  studies  of  flow  around  the  DTM  B  5415 


Figure  1.  Grey  scale  image  of  the  DTMB  model  5415  geometry 


The  DTMB  model  5415  was  first  studied  extensively  using  CFD  techniques  when  it  was 
chosen  to  beoneof  three  test  cases  in  the  Gothenburg  2000  workshop  for  computational  fluid 
dynamics  applied  to  ship  flows.  Participants  were  asked  to  focus  in  particular  on  the  total 
resistance,  the  wave  profile  along  the  hull  and  at  selected  locations  close  to  the  ship,  the 
overal  I  w  ave  pattern  i  n  the  near  fi  el  d  and  the  mean  f  I  ow  vel  oci  ti  es  near  the  stern ,  parti  cu  I  arl  y 
at  the  propeller  plane.  Experimental  data  for  this  model  has  been  obtained  by  three 
laboratories;  DTMB,  1 1 FH  R,  and  INSEAN  (Instituto  Nazionaleper  Studi  ed  Esperienzedi 
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A  rchi  tettura  N  aval  e).  A  detai  I  ed  compari  son  of  thi  s  data  i  s  provi  ded  i  n  the  paper  by  Stern  et 
al  .[23].  Figure  1  shows  a  gray-scale  image  of  the  model  geometry  and  thetest  conditions  for 
the  simulation  are  given  in  Table! 

Seven  CFD  codeswereused  to  model  theDTM  B  5415attheGothenburg  2000 workshop;  CFX, 
COMET,  FINFLO,  ICARE,  CFDSHIP,  MGSHIP  and  UNCLE.  The  main  conclusions  drawn 
from  theanalyses  of  the  seven  si  mul  ati  ons  and  a  compari  son  with  theexperi  mental  data  were 
as  follows: 

1.  Thecomputed  total  drag  coefficient  (Ct)  differed  from  the  experimentally  measured  value 

by  a  maxi  mum  of  ±  7%,  and  theaveragevalueovertheseven  codes  was  al  most  exactly  the 
same  as  the  data  (0.5%  higher). 

2.  M  ost  of  the  codes  had  difficulty  accuratel  y  si  mu  I  ati  ng  the  wave  prof  i  I  es  i  n  the  near  fi  el  d . 
Themain  problem  was  theunder  prediction  of  the  peak  in  thebow  wavedueto  excessive 
numerical  damping. 

3.  The  flow  in  the  propeller  plane  was  reasonably  well  predicted  by  half  of  the  codes.  No 
vortex  could  be  detected  in  the  cross-flow  vector  plots,  but  detailed  vorticity  plots  were 
able  to  capture  the  vortex  shed  from  the  bulbous  bow. 


Table  1.  Test  Conditions  for  scale  model  at  20  knots 


Scale  Ratio 

24.832 

Length  (L) 

5.72  m 

Draft  (T) 

0.248  m 

Wet  Surface  Area  (S) 

4.861  m2 

AdvanceVelocity  (Uo) 

2.0637  m/  s 

Froude Number  (Fn) 

0.2755  m /  s 

Reynolds  N  umber  (Re) 

1.26  x  107 

Sinkageat  FP 

-0.0028L 

Sinkageat  AP 

-0.0009L 

With  regard  to  the  flow  into  the  propeller  plane  for  the  DTBMB  5415,  Gorski  [24]  has 
commented  "The  dominant  flow  feature  is  a  vortex  created  at  the  sonar  dome  that  flows  downstream 
to  the  propeller  plane.  In  general...  RANS  computations  with  two-equation  models  can  predict  this 
vortical  flow  well.  The  exact  position  and  strength  of  the  vortex,  which  depends  on  grid  resolution  and 
turbulence  modelling,  are  not  necessarily  obtained,  but  enough  detail  of  the  average  inflow  and  wake 
deficits  for  propeller  design  are  obtained.  At  full  scale  the  vortex  will  be  closer  to  the  hull  and  may  or 
may  not  impact  the  propeller  significantly." 

Subsequent  to  the  Gothenburg  2000  workshop  the  DTM  B  model  5415  has  been  used  as  a  test 
case  for  a  number  of  other  CFD  codesand  interface  model  ling  techniques.  Rheeand  Skinner 
[19]  used  Fluent  version  6. 2  to  si  mu  I  ate  flow  around  theDTM  B  5415  under  thetest  conditions 
shown  in  Table!  They  used  a  half  domain  hexahedral  cell  mesh  containing  approximately 
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1.8  mi  1 1  ion  cel  Is  and  the  HRIC  (High  Resolution  InterfaceCapturing)  VOF  formulation  within 
Fluent  to  capture  the  air/  water  interface.  Simulations  were  performed  using  two  different 
turbul  ence  model  s;  the  shear-stress-transport  k-co  model  of  M  enter  and  the  Reynol  ds  stress 
transport  model.  The  total  resistance  coefficient  was  under  predicted  by  9.4%  but  the  wave 
profiles  along  the  hull  surface  and  in  the  near  field  agreed  very  well  with  the  measured  data. 
Peak  values  were  under  predicted  away  from  the  hull  however  and  this  was  attributed  to 
numerical  diffusion  enhanced  by  stretched  cells  in  the  area  outside  of  the  imminent 
neighbourhood  of  the  hull.  The  velocity  components  in  the  propeller  plane  showed  good 
agreement  with  the  experimental  values. 

Wood  et  al.  [25]  used  the  finite  volume  commercial  solver  CFX  to  perform  a  similar 
calculation  for  the  DTM  B  5415.  The  test  conditions  were  again  the  same  as  those  shown  in 
Table!  They  used  the  VOF  model  implemented  in  CFX  to  model  the  air/  water  interface  and 
the  shear-stress-transport  k-co  model  fortheturbulentflow.  Simulations  were  performed  using 
both  structured  hexahedral  meshes  on  a  half  domain  mesh  or  unstructured  tetrahedral 
meshes  with  prism  boundary  layersonthefree  surface  in  an  attempt  to  obtain  accurate  wave 
profile  simulations.  Results  calculated  on  the  unstructured  meshes  however  showed  very 
poor  agreement  with  the  experimental  results.  On  the  structured  hexahedral  mesh  they 
calculated  the  drag  coefficient  to  have  the  val  ue  4.368  x Id3,  which  is  3.2%  higher  than  the 
experimental  valueof  4.23xlOr3  [3].  Waveprofiles  along  thehull  and  in  the  near  field  were 
also  obtained  and  in  general  showed  good  agreement  with  the  experimental  results,  although 
their  calculated  profiles  also  showed  a  small  loss  in  wave  amplitude  in  the  peaks  of  the 
p  rof  i  I  es,  w  h  i  ch  th  ey  attr  i  buted  as  bei  n  g  d  u  e  to  n  u  mer  i  caldiffusion.Axialvelocityplotsinthe 
propeller  plane  also  showed  good  agreement  with  the  experimental  vales. 

D  i  M  asci  o  et  al .  [  15]  have  used  the  DTM  B  5415  to  val  i  datethei  r  si  ngl  e-phase  I  evel  set  method 
coupled  with  a  standard,  in-house,  CFD  RANS  solver.  Their  code  uses  the  finite  volume 
techniquewith  pressure  and  velocity  co-  located  atthecell  centre  and  a  second  order  Godunov 
type  algorithm  to  discretize  the  RANS  equations.  The  Spalart-Allmaras  one-equation 
turbul  ence  mod  el  was  used  fortheturbulentflow.  Thephysi  cal  domain  was  meshed  usinga 
structured  hexahedral  grid  with  approximately  2.5million  cells.  For  thetest  conditions  given 
inTablelthey  calculated  atotal  drag  coefficient  of  4.39xlOr3,  which  is  3.8%  higher  than  the 
experimental  value.  Their  computed  wave  shape  along  the  hull  shows  excellent  agreement 
with  the  experimental  data,  although  there  is  a  slight  reduction  in  peak  value  at  the  bow 
compared  to  the  experi  mental  vale,  which  Di  Mascio  et  al.  propose  can  be  explained  by  the 
way  in  which  the  photographic  and  digitizing  system  was  used  to  obtain  the  experi  mental 
results. 


4.  The  Fluent  Code 

Fluent  version  6.3  is  a  cell  centred  finite  volume  general  purpose  CFD  computer  code 
developed  by  Fluent,  Inc.  and  now  marketed  by  ANSYS,  Inc.  an  engineering  simulation 
software  provider  which  acquired  Fluent,  Inc.  in  2006.  The  Hydrodynamics  Research  Group 
within  M  PD  has  been  using  this  code  for  several  years  to  simulate  a  number  of  problems 
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involving  submarine  related  flows  [7, 8, 9, 10],  This  report  describes  our  first  application  of  the 
code  to  surface  ship  flows. 

A  disti  nctdifference  between  the  si  mulationsdescri  bed  hereand  thosereported  earlieristhe 
presence  of  the  air/  sea  interface  in  the  calculation.  Asexplained  in  Section  2,  there  are  two 
main  approaches  to  maintaining  a  distinct  interface  in  a  simulation,  either  surface  fitting 
approaches  or  surface  capturing  approaches.  Fluent  has  implemented  the  surface  capturing 
approach  by  using  theVOF  scheme  for  general  multiphase  flow  modelling.  This  involves 
defining  a  volume  fraction  function  for  each  of  the  fluids  throughout  thedomain  and  then 
converting  the  volume  fraction  of  each  fluid  with  the  average  fluid  flow.  The  interface 
betw  een  the  tw  o  f I  u  i  d  s  i  s  then  reconstru  cted  from  the  vol  u  mef  radii  on  f u  ncti  on  for  each  of  the 
fluids  in  the  vicinity  of  the  i  nterface. 

Fluent  6.3  offers  a  choice  of  several  different  methods  for  interface  reconstruction.  The 
appropriate  scheme  to  use  depends  on  whether  the  simulation  is  steady-state  or  time- 
dependent  and  whether  an  implicitorexplicitti  me  discretization  schemeisused.Theresults 
described  in  this  report  were  calculated  from  implicit  steady-state  simulations  and  a  H  igh 
Resolution  I  nterface  Capturing  (FI  RIC)  schemewasused  to  ensure  a  sharp  interface  betw  een 
thetwo  fluids.  This  is  a  modification  of  the  FI  Rl  C  schemedescri  bed  by  Muzaferijaetal.  [26],  It 
consists  of  anon-linear  blend  of  upwind  and  downwind  discretizationssothatthecomputed 
fluxes  ofthevolumefractionsdonotunderflow  or  overflow  thecells.  Standard  interpolation 
schemes  are  used  to  obtai  n  thefacefl  uxes  w  henever  a  cel  I  i  s  compl  etely  fi  1 1  ed  wi  th  onefl  ui  d . 

Simulating  the  motion  of  a  surface  ship  not  only  requires  maintaining  a  sharp  interface 
between  the  sea  and  air  phases  but  also  requires  specification  of  appropriate  boundary 
conditions  at  the  inlet  and  outlet  of  thedomain  for  each  of  thefluids  in  the  simulation.  One 
way  to  do  this  is  to  write  a  User-Defined  Function  to  specify  the  total  pressure  and  VOF 
profile  for  each  of  thefluids  at  both  the  inlet  and  outlet.  The  other  option  is  to  use  the  new 
Open  Channel  Boundary  Condition  implemented  in  Fluent  6.3.  Using  this  method  the 
i  nf  I  ow/  outf  I  ow  cond  i  ti  on  i  s  sped  f  i  ed  vi  a  the  i  nf  I  ow  vel  oci  ty  and  thef  ree  su  rf  ace  I  ev el  and  the 
pressureand  VOF  profiles  are  automatically  calculated.  This  simplifies  theinitial  set-up  of  the 
simulation  and  was  the  boundary  condition  used  in  the  simulations  described  here. 

Fluent  6.3  offers  a  variety  of  different  RAN  S  turbulence  models.  These  include  the  Spalart- 
Allmaras  model,  thefc-s  model,  and  thefc-co  model,  all  of  which  are  based  on  theBoussinesq 
approach,  which  assumes  that  the  calculated  turbulent  viscosity  coefficient  is  isotropic. 
FI  ow  ever  for  compl  ex  f  I  ow  s  i  nvol  vi  ng  streaml  i  ne  cu  rvatu  re,  sw  i  rl ,  rotati  on,  and  rapi  d  changes 
i  n  strai  n  rate  thi  s  i  s  not  the  case.  For  these  fl  ows  FI  uent  provi  des  the  Reynol  ds  Stress  M  odel 
(RSM ),  which  solves  the  RA  NS  equations  by  solving  additional  transport  equations  for  each  of 
the  individual  Reynolds  stresses.  This  means  thatfiveadditional  transport  equations  must  be 
solved  in  three-dimensional  flows  and  the  simulation  timeincreases  substantially.  Si  nee  the 
si  mul  ati  ons  descri  bed  i  n  thi s  report  i  nvol  ve onl y  strai  ght-ahead  moti on  at  mi  ni  mal  angl es  of 
attack  weexpecttheflow  to  remain  largely  attached  tothehull,  hencewehaveemployed  the 
standard  k-s  turbulence  mod  el  and  standard  wall  fundi  ons  for  these  si  mu  I  ati  ons.  Rheeand 
Skinner  [19]  have  previously  used  both  M  enter's  shear-stress  transport  model  and  an  RSM 
model  to  simulate  the  same  test  case  as  that  considered  here  and  have  shown  that  the 
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boundary  layer  remains  attached  to  the  hull  for  the  majority  of  the  flow  and  issimilar  to  an 
equilibrium  boundary  layer  flow. 


5.  Simulation  Results  with  a  Fully  Hexahedral  Mesh 

A  fully  structured  hexahedral  mesh  containing  approximately  3.8  million  cells  was 
constructed  using  the  Fluent  preprocessor,  Gambit.  As  the  flow  has  a  plane  of  symmetry 
about  the  centre  plane  a  half  domain  grid  was  used.  Thehull  surface  was  first  meshed  using 
quadrilateral  elements.  A  volume  forward  of  the  bow  was  then  created  by  projecting  the 
outli  ne  of  the  bow  surface  towards  the  inlet.  The  mesh  for  this  vol  u  me  was  created  from  the 
projection  of  the  surface  mesh  on  the  matching  portion  of  the  bow.  Projection  of  the  stern 
geometry  and  its  surface  mesh  towards  the  outlet  created  the  corresponding  volume  and 
mesh  for  the  stern.  The  i  ni  ti  al  secti  on  of  the  bow  projecti  on  was  angl  ed  to  mi  ni  mi  se  the  skew 
of  the  el  ements.  Th  i  s  w  as  u  n  necessary  for  the  stern  geometry  as  the  ori  entati  on  of  the  stern  i  s 
approximately  parallel  to  that  of  the  outlet  and  thus  there  is  no  tendency  for  highly  skewed 
el  ements  to  occur  in  this  region. 

The  deck  of  the  vessel  and  theupper  surface  of  the  volumes  projected  from  the  bow  and  stern 
were  i  n  turn  projected  to  the  desi  red  hei  ght  for  the  top  of  thefl  ui  d  domai  n,  thus  creati  ng  the 
vol  u  mes  above  these  su  rf  aces.  The  mesh  for  these  vol  u  mes  w  as  created  from  the  p  roj  ecti  on  of 
these  surface  meshes.  A  "C"  volume  and  grid  in  thevertical  plane  was  wrapped  around  the 
hull  and  projected  volumes.  The  grading  of  the  elements  around  the  bow  was  adjusted  to 
minimise  the  skew  ness  around  thesonar  dome.  Carewastakentoensureasuitabledensity  of 
el  ements  at  the  waterline,  bow,  stern,  and  wakeregion.  Details  of  the  surface  mesh  over  the 
bow  are  shown  in  Figure  2,  and  the  surface  mesh  over  the  entire  hull  isshown  in  Figure  3. 

The  experimental  conditions  for  the  simulation  were  the  same  as  those  in  the  Gothenburg 
2000  workshop  and  are  shown  in  Table  1.  The  attitude  of  the  model  with  respect  to  the 
coordinate  axes  was  set  according  to  the  experimentally  measured  sinkageand  trim  values 
before  meshi  ng  commenced .  The  ori  gi  n  of  the  coordi  nate  system  was  I  ocated  at  the  mi  dshi  p 
intersection  of  the  calm  water  free-surface  and  the  centre  plane.  The  open  channel  boundary 
condition  was  used  to  specify  the  inlet  and  outlet  boundary  condition,  a  symmetry  planewas 
used  along  the  centre  plane,  and  the  remaining  boundary  surfaces  along  the  exterior  of  the 
domain  were  set  to  slip  wall  conditions.  The  longitudinal  extent  of  the  mesh  was-  2Lpp<x< 
5Lpp  and  the  radial  extent  was  0.0  <  r  <  2.5Lpp.  (Lpp  denotes  the  length  between 
perpendiculars). 
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Figure  2.  Detailed  view  of  the  surface  mesh  on  the  bidbous  bow 

An  implicit  steady-state  cel  I  based  solution  procedure  was  used  to  solve  the  equations.  The 
SIMPLE  algorithm  was  used  for  the  pressure-velocity  coupling,  the  PRESTO  scheme  for  the 
pressure  interpolation,  the  2nd  order  upwind  scheme  for  the  solution  of  the  momentum 
equati  ons  and  the  mod  if  i  ed  H  Rl  C  scheme  for  the  sol  uti  on  of  the  vol  u  me  f  racti  on  equati  ons. 
The  relaxation  factors  were  typically  set  to  0.2.  The  standard  k- s  turbulence  model  with 
equi  I  i  bri  um  wal  I  functi  ons  was  used  to  si  mu  I  ate  turbul  ent  fl  ow .  The  y+  val  ues  for  the  wal  I 
adjacent  cells  over  the  hull  were  in  the  range  80.0  <y+<  100.0,  which  is  comfortably  within 
the  guidelines  for  the  use  of  the  wal  I  function  approach. 

Initial  convergence  of  the  simulation  was  found  to  be  considerably  enhanced  by  paying 
careful  attention  to  the  initialization  procedure.  Theprimary  phase  should  always  be  set  to  the 
lower  density  fluid,  which  in  this  case  is  air.  The  specified  operating  density  should  be  set  to 
that  of  the  pri  mary  phase  and  the  reference  pressure  I  ocati  on  shou  I  d  be  set  to  a  regi  on  w  hi  ch 
willalwayscontaintheprimaryph  ase.  1 1  w  as  al  so  f  ou  n  d  h  el  pf  u  I  to  i  n  i  ti  al  i  ze  the  enti  re  w  ater 
domain  with  the  correct  hydrostatic  pressure  profile  and  to  initialize  both  the  water  and  air 
domai  n  to  the  same  vel  ocity. 
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Figure  3.  Surface  mesh  over  the  entire  hull 

I  ni  ti  al  si  mu  I  ati  ons  u  si  ng  thi  s  mesh  w  ere  d  i  scou  ragi  ng.  A  fter  13,000  i  terati  ons  the  w  ater  su  rf  ace 
showed  no  signs  of  developing  the  classical  Kelvin  wake  pattern,  even  though  all  residuals 
had  decreased  by  several  orders  of  magnitude  at  this  stage.  M  ore  importantly,  the  water 
surface  contours  indicated  the  presence  of  significant  computational  errors  occurring  at  the 
inlet  and  at  thefar  field  boundary.  Thesourceof  these  errors  was  not  definitively  identified 
but  it  was  thought  that  they  originated  either  from  an  incompatibility  of  theOpen  Channel 
boundary  condition  with  thecurved  outer  boundary  surface,  or  from  the  irregular  natureof 
the  surface  mesh  over  the  inlet  boundary  due  to  the  meshing  constraints  imposed  by  the 
projection  of  the  bulbous  shape  of  the  lower  bow.  The  nature  of  these  errors  is  indicated  in 
Figure4,  which  shows  a  contour  plot  of  the  vertical  height  of  the  water  surface,  which  isthe 
iso-surface  of  all  cellswherethevolume  fraction  ofwater  has  thevalue0.5.  Thecurved  nature 
of  thefar  field  geometry  is  also  evident  in  Figure4,  which  shows  the  surface  mesh  over  the 
upstream  inlet  boundary.  The  physical  extent  of  thedomain  isevidentfromtheoutlineofthe 
hull. 

As  a  fi  rst  attempt  to  obtai  n  a  real  istic  Kel  vi  n  wake  pattern  it  was  deci  ded  to  i  ncreasethe  mesh 
resolution  in  the  vicinity  of  the  water  surface.  To  do  this  the  "adaption"  feature  of  the 
softwarewas  used,  which  automatically  refi  nes  the  mesh  in  a  given  region  by  splitting  each 
hexahedral  cell  into  eight  smaller  hexahedral  cells.  Whilst  being  easy  to  implement,  this 
procedure  results  in  a  significant  increase  in  the  total  number  of  cells  on  the  mesh.  An 
adapti  on  w  as  then  performed  i  n  the  regi  on  of  the  mesh  20  cm  above  and  bel  ow  the  water  I  i  ne 
and  thi  s  took  the  total  nu  mber  of  cel  I  s  to  approxi  matel  y  9.5  mi  1 1  i  on .  The  si  mu  I  ati  on  was  then 
allowed  to  run  out  to  35,000  iterations  but  the  results  were  again  discouraging  as  no 
i  improvement  was  seen  i  n  the  shape  of  the  water  surface. 
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Figure  4.  Water  surface  elevation  and  nature  of  curved  outer  boundary 


A  real  istic  Kelvin  wake  pattern  was  fi  nal  I  y  obtai  ned  by  add  ressi  ng  the  probl  ems  descri  bed  i  n 
the  previ  ous  paragraph.  Two  changes  were  made:  (a)  the  curved  outer  boundary  surface  was 
removed  and  replaced  by  a  rectangular  box  geometry.  The  resulting  volume  was  then  meshed 
using  regular  hexahedral  cells,  (b)  Additional  hexahedral  cells  were  added  upstream  of  the 
exi  sti  ng  mesh  to  provi  de  a  uniform  surface  mesh  for  the  i  mpositi  on  of  the  i  nfl  ow  boundary 
conditions.  A  non-conformal  interface  was  used  between  thefaceofthisadditional  mesh  and 
the  irregular  surface  mesh  at  the  original  mesh  boundary.  These  changes  increased  the 
number  of  cel  Is  in  the  mesh  to  approximately  5.6  mil  lion.  The  new  mesh  was  then  run  out  to 
8, OCX)  iterations  until  the  calculation  converged.  It  should  be  noted  that  our  criterion  for 
convergence  is  based  on  a  combination  of  monitored  metrics  and  a  reduction  in  calculated 
residuals.  In  the  simulations  reported  herewemonitored  theshapeofthefree  surface  as  well 
as  the  reduction  in  residuals.  We  considered  the  calculation  to  have  converged  when  an 
increase  in  the  number  of  iterations  produced  no  further  change  to  the  shape  of  the  free 
surface  and  left  the  magnitude  of  the  residuals  unchanged.  The  resulting  wave  pattern  is 
shown  in  Figure  5,  and  Figure  6  shows  a  more  detailed  view  of  the  wave  pattern  in  the 
vicinity  of  the  hull. 
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Figure  5.  Wake  pattern  obtained  after  8,000  iterations  using  a  rectangular  domain 


Figure  6.  Detailed  view  of  the  wave  pattern  in  the  vicinity  of  the  hull 
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The  total  resistance  coefficient  calculated  from  this  simulation  isthesum  of  thetotal  viscous 
force  plus  the  total  pressure  force  on  the  hull  and  has  the  value  4.42  xlO'3,  which  is  4.5% 
higher  than  the  experimental  value  of  4.23  xlO-3  [3].  The  friction  coefficient  has  the  value 
3.07  xlO'3,  which  is  only  1.3%  higher  than  the  average  result  from  the  seven  simulations 
reported  at  the  Gothenburg  2000  workshop  [3],  while  the  pressure  coefficient,  which  hasa 
val  ue  of  1.36  x  10-3 ,  i  s  12.4%  hi  gher  than  the  average  val  ue  from  the  same  seven  si  mu  I  ati  ons. 
Rheeand  Skinner  [19],  who  more  recently  performed  a  si  mu  I  ati  on  of  theDTM  B  hull  using  the 
Fluent  code  and  the  more  sophisticated  RSM  turbulence  model,  found  a  total  resistance 
coefficient  of  3.83  xlO'3,  which  isan  under  prediction  of  9.4%  from  the  experi  mental  value. 

Figure  7  shows  the  simulated  wave  profile  along  the  hull  surface  compared  with  the 
experimental  data  from  the  Gothenburg  test  set  [27].  Overall  there  is  good  qualitative 
agreement  between  the  simulated  and  experimental  profiles.  The  phase  and  amplitudeofthe 
primary  bow  wave  are  slightly  in  error  but  this  is  a  common  problem  in  this  type  of 
simulation.  Rautaheimo  and  Salminen  [28]  under  predicted  the  peak  amplitude  by  25%  in 
their  simulation  of  the  DTMB  5415  using  the  FINFLOW  code.  The  only  significant 
disagreement  with  the  experi  mental  data  isthestrong  dip  in  thewater  surface  immediately 
foil  owing  the  bow  wave.  This  is  not  seen  inthesi  mul  ati  ons  of  Rautaheimo  and  Salminen  [28], 
Rhee  and  Ski  nner  [  19]  or  C  hen  et  al .  [  18],  but  i  s  noti  ceabl  e  i  n  the  si  mul  ati  on  of  Li  [  14]  usi  ng  the 
FI  NFLOW-RA  NS  SHIP  solver. 

Better  agreement  between  thesimulated  and  experimental  data  can  be  seen  in  Figure8,  which 
shows  the  wave  profi  I  e  al  ong  the  I  i  ne  at  y/  Lpp  =  0.172.  There  i  s  excel  I  ent  agreement  i  n  the 
shape  of  the  profile,  although  the  phase  difference  between  thesimulated  and  experimental 
results  slowly  increases  along  the  length  of  the  cut. 

The  ability  to  accurately  simulate  the  wake  flow  in  the  vicinity  of  the  propeller  plane  is 
i  important  for  the  pred  i  cti  on  of  propel  I  er  I  oad  i  ngs  and  Fi  gu  re  9  shows  a  compari  son  betw  een 
thesimulated  and  experimental  velocity  components  on  thepropeller  planeatz/  Lpp  =-0.02. 
There  is  excellent  agreement  for  they  (v/  U)  and  z  (w/  U)  components  whilethe  maximum 
discrepancy  for  the  axial  component  (x/  U)  is  around  10%. 

To  check  any  grid  dependency  of  these  results  weadapted  the  grid  in  the  region  20  cm  above 
and  below  thewater  line,  resulting  in  a  mesh  with  approximately  11.5  million  hexahedral 
cel  I  s.  Thesi  mul  ati  on  was  then  conti  nued  out  to  17,000  iterati  ons  unti  I  the  resi  dual  s  had  agai  n 
decreased  to  acceptable  levels.  Visual  inspection  of  thewater  surface  showed  effectively  no 
change  to  the  wave  pattern  and  detailed  comparison  with  theindividual  wave  profiles  and 
the  results  calculated  on  the  coarser  mesh  showed  that  there  was  very  little  change  to  the 
simulated  results. 
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Figure  7.  Comparison  of  simulation  with  experimental  results:  wave  profile  along  the  hull  for  fully 
structured  hexahedral  mesh 


Figure  8.  Comparison  of  simulation  with  experimental  results:  wave  profile  along  the  line  at  yA-PP  = 
0.172  for  fully  structured  hexahedral  mesh 

Plots  for  the  wave  profile  along  the  hull  and  along  the  line  at  y/  Lpp  =0.172  are  shown  in 
Figures  10  and  11  respectively.  The  finer  mesh  has  resulted  in  a  slight  improvement  in  the 
phase  of  the  peak  bow  wave,  but  the  magnitude  of  the  peak  has  moved  si  i  ghtly  further  from 
the  experimental  result.  A  similar  result  has  occurred  for  the  wave  profile  along  the  line  at 
y/  Lpp  =0.172,  where  thefiner  mesh  has  resulted  in  an  increased  height  for  the  primary  wave, 
but  this  has  moved  the  calculated  value  further  from  the  experimental  result. 

Thetotal  resistance  coefficient  on  theunadapted  mesh  has  the  value  4.42  xlO'3,  which  is4.5% 
higher  than  the  experimental  value.  On  the  adapted  mesh  this  changes  to  4.39  xlO'3,  which 
agrees  slightly  better  with  the  experimental  value  of  4.23xlOr3[3]  and  isonly  3.8%in  error. 
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Figure  9.  Comparison  of  simulation  with  experimental  results:  velocity  components  on  propeller  plane 
at  z/Lpp  =  -0.02  for  fully  structured  hexahedral  mesh 

Fi  gure  12  shows  the  effect  of  the  gri  d  refi  nement  on  the  vel  oci  ty  components  i  n  the  propel  I  er 
plane.  For  they  and  z  components  there  is  effectively  no  change  to  the  simulated  results, 
while  for  the  axial  component  the  increased  grid  resolution  has  resulted  in  a  very  slight 
change  i  n  the  resu  1 1  near  y/  L  p  p  =  0. 015  which  bri  ngs  it  just  si  i  ghtly  cl  oser  to  the  experi  mental 
result. 


Figure  10.  Effect  of  grid  refinement:  wave  profile  along  the  hull  for  original  hexahedral  mesh  and 
refined  hexahedral  mesh 
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Figure  11.  Effect  of  grid  refinement:  wave  profile  along  the  line  at  y/Lpp  =  0.172  for  original  hexahedral 
mesh  and  refined  hexahedral  mesh 


Figure  12.  Effect  of  grid  refinement:  velocity  components  on  propeller  plane  at  z/Lpp  =  -0.02  for 
original  hexahedral  mesh  and  refined  hexahedral  mesh 
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6.  Simulation  Results  with  a  Tetrahedral/Prism  Mesh 

One  of  the  advantages  of  the  FI  uent  code  for  si  mu  I  ati  ons  of  naval  vessel  s  is  its  ability  to  solve 
the  RA  NS  equations  on  either  a  structured  hexahedral  mesh  or  an  unstructured  hybrid  mesh 
consisting  of  tetrahedral  cellsand  prism  cells.  It  is  generally  accepted  that  hexahedral  meshes 
providemoreaccuratesi  mulation  results[25],  buttheonedisadvantageofthisapproach  isthe 
time  required  to  construct  the  mesh.  The  hexahedral  mesh  used  in  the  previous  section  took 
approximately  five  days  to  construct,  while  the  hybrid  mesh  used  in  this  section  was 
constructed  in  oneday.  Whilst  the  unappended  DTMB  5415  is  a  relatively  simpleshapeto 
mesh  using  either  a  hexahedral  or  hybrid  scheme  this  will  not  necessarily  be  the  case  when 
simulations  on  fully  appended  hulls  are  required.  In  that  case,  especially  when  using  the 
Gambit  software  package  to  create  the  mesh,  the  use  of  an  unstructured  hybrid  mesh  will 
resultin  considerabletime  savings.  In  this  section  we  have  used  an  unstructured  hybrid  mesh 
to  gain  some  experience  with  this  technique  in  anticipation  of  the  need  to  consider  more 
complicated  geometries  in  future  naval  applications. 

A  hal  f  d  omai  n  gri  d  w  as  agai  n  used  but  the  si  ze  of  the  d  omai  n  w  as  smal  I  er  than  that  u  sed  w  i  th 
the  hexahedral  mesh,  i  n  this  case  extending  a  distance  2  Lpp  aft  of  thehul  I,  one  Lpp  upstream 
from  the  bow,  1.5  Lpp  out  from  the  symmetry  plane,  and  the  total  vertical  depth  of  the 
domain  was  one  Lpp.  The  undisturbed  sea  surface  was  located  on  thez  =0  plane.  A  surface 
mesh  for  the  hull  was  constructed  from  triangular  elements  with  side  lengths  of 
approximately  2cm.  A  Gambit  curvature  size  fundi  on  was  used  to  ensuresmooth  growth  of 
the  mesh  si  ze  i  n  hi  ghly  curved  regi  ons  of  the  hull.  The  maxi  mum  angl  e  between  the  normal  s 
to  the  surface  was  set  to  20°.Thehull  top,  undisturbed  sea  surface  and  symmetry  planewere 
meshed  with  triangular  elements  having  sizes  of  5cm,  10cm  and  50cm  respectively  with 
growth  rates  of  approximately  1.15.  Prism  layers  were  then  grown  from  the  hull  surface, the 
symmetry  pi  ane  and  the  pi  ane  represent  ng  the  u  nd  i  sturbed  sea  su  rface.  The  fi  rst  cel  I  hei  ght 
for  each  of  these  surfaces  was  set  to  5.72mm.  This  provides  good  resolution  for  waves  on  the 
sea  su  rface  but  i  s  far  too  I  arge  to  ach  i  eve  the  opti  mu  m  y + val  u  e  of  100  on  the  shi  p  su  rface.  W  e 
therefore  expect  that  our  simulated  friction  drag  will  not  be  very  accurate,  but  the  wave 
profiles  should  bewell  resolved.  A  tetrahedral  volume  mesh  was  then  grown  from  each  of 
these  surfaces,  again  using  appropriate  size  fundi  ons  to  ensure  a  uniform  rate  of  change  in 
cell  size,  and  the  resulting  mesh  had  1.8  million  cells. 

The  simulation  was  again  run  using  an  implicit  steady-state  cell  based  solution  procedure 
with  the  SIMPLE  algorithm  for  the  pressure-velocity  coupling,  the  PRESTO  scheme  for  the 
pressure  interpolation,  the  2nd  order  upwind  scheme  for  the  solution  of  the  momentum 
equations  and  the  modified  HRIC  scheme  for  the  solution  of  the  volume  fradi  on  equations. 
The  standard  k-e  turbulence  model  with  equilibrium  wall  fundi  ons  was  again  used  to 
si  mu  I  ate  tu  rbu  I  ent  f  I  ow .  The  si  mu  I  ati  on  w  as  ru  n  for  25,000  i  terati  ons  and  the  resu  I  ti  ng  w  ave 
pattern  is  shown  in  Figure  13. 

The  total  resistance  coefficient  from  this  simulation  has  the  value  3.40  xHT3,  which  is 
approxi  mately  20%  lower  than  the  experi  mental  value  of  4.23  xlO- 3  [3].  This  is  quite  a  good 
result,  consideringthatthey+valuearound  the  hull  surface  is  approxi  mately  1.0  xlO4,  which 
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i  s  two  orders  of  magnitude  greater  than  the  opti  mum  val  uefor  the  use  of  wal  I  functi  ons.  The 
resu  1 1  i  s  somew  hat  fortu  i  tou  s  how  ever  as  thef  ri  cti  on  coeff  i  ci  ent  has  the  val  ue  1.6  x  lCf3,  w  hi  ch 
is  approximately  40%  lower  than  the  average  result  from  the  simulations  reported  at  the 
Gothenburg  2000  workshop  [3].  This  lower  frictional  drag  is  roughly  compensated  by  a  higher 
pressuredrag,  si  ncethe  pressure  coefficient  has  a  value  of  1.80xl0-3,  which  is  approximately 
48%  hi  gher  than  the  average  val  ue  from  the  si  mul  ati  ons  at  the  workshop. 


z-cocadmaie 


0.025 

0.02 

0.015 


0.01 

0.005 

0 

•0.005 

•0.01 


Figure  13.  Wave  pattern  obtained  after  25,000  iterations  on  a  tetrahedral  mesh. 

Figure  14  shows  the  simulated  wave  profile  along  the  hull  surface  compared  with  the 
experimental  data  from  theGothenburg  test  set  [27]  and  the  previous  si  mul  ati  on  results  from 
the  (unadapted)  hexahedral  mesh.  Once  again  there  is  overall  good  qualitative  agreement 
betw  een  th e  si  mu  I  ated  an d  exp er i  mental  p  rof  i  I  es,  w  i  th  th e  d  eg ree  of  er ror  f  rom  th e  resu  I  ts  on 
the  hybrid  mesh  being  similar  to  those  from  the  hexahedral  mesh.  The  phase  of  the  primary 
bow  wave  is  accurately  predicted  buttheamplitudeistoolow.Thisisundoubtedly  duetothe 
lower  resolution  ofthehybrid  mesh  compared  to  the  hexahedral  mesh  in  the  vicinity  of  the 
sea/  ai  r  i  nterface.  The  strong  d  i  p  i  n  the  water  surface  i  mmed  i  atel  y  fol  I  ow  i  ng  the  bow  wave  i  s 
sti  1 1  present  but  i  s  not  as  pronou  need ,  although  this  coul  d  agai  n  be  d  ueto  the  I  ower  resol  uti  on 
ofthehybrid  mesh. 

Figure  15  showsthewaveprofilealongthelineaty/  Lpp  =0.172  and  it  can  be  seen  that  there 
is  quite  good  agreement  between  the  profiles  calculated  on  the  hybrid  mesh  and  on  the 
hexahedral  mesh.  The  major  difference  between  the  two  results  again  being  evident  in  the 
calculated  amplitudes,  where  the  relative  coarseness  of  the  hybrid  mesh  has  resulted  in 
significantly  lower  peak  amplitudes. 

The  calculated  velocity  components  on  the  propeller  plane  at  z/  Lpp  =-0.02  are  shown  in 
Figure  16.  The  agreement  between  the  y  component,  (v/  U ),  computed  on  the  hybrid  mesh 
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and  on  the  hexahedral  mesh  is  excellent,  and  both  results  are  in  good  agreement  with  the 
experimental  result.  Forthez  component,  (w/  U),  there  is  a  slight  disagreement  between  the 
results  simulated  onthehybrid  and  thehexahedral  meshes  quiteclose  to  thesymmetry  plane 
but  the  results  soon  becomeindistinguishableasy/  Lpp  increases.  The  axial  component(x/  U) 
displays  quitea  different  shapefrom  both  the  experimental  data  and  the  simulated  result  on 
the  (unadapted)  hexahedral  mesh  for  y/  Lpp  <  0.04,  but  good  agreement  over  the  rest  of  the 
range. 


Figure  14.  Comparison  of  simulation  results  computed  on  hexahedral  and  tetrahedral  meshes:  wave 
profile  along  the  hull 
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Figure  15.  Comparison  of  simulation  results  computed  on  both  hexahedral  and  tetrahedral  meshes: 
wave  profile  along  the  line  at  y/Lpp  =  0.172 


Figure  16.  Comparison  of  simulation  results  computed  on  both  hexahedral  and  tetrahedral  meshes: 
velocity  components  on  propeller  plane  at  z/Lpp  =  -0.02 

To  check  the  effect  of  the  gri  d  on  these  resu  I  ts  we  agai  n  adapted  the  mesh  i  n  the  regi  on  20  cm 
above  and  below  the  water  line.  The  adapted  mesh  had  3.48  million  mixed  cells  and  the 
simulation  was  then  continued  outto  38,500  iterations  until  the  calculation  converged.  Figure 
17  shows  the  effect  of  the  gri  d  refi  nement  on  the  wave  profi  I  e  al  ong  the  hull.  Thefi  ner  mesh 
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has  resulted  in  an  increase  in  the  amplitude  of  the  bow  wave  so  that  it  now  agrees  with  the 
experi  mental  resu  1 1.  A  part  from  the  i  ncrease  i  n  peak  ampl  i  tude  how  ever  thi  s  i  s  essenti  al  I  y  no 
difference  in  the  profiles  calculated  on  thetwo  meshes. 

A  similar  result  has  occurred  for  the  wave  profile  along  thelineaty/  Lpp  =0.172,  shown  in 
Figure  18.  Thefiner  mesh  has  resulted  in  an  increased  height  for  the  primary  and  secondary 
waves  and  this  has  moved  thecalculated  result  closer  to  the  experi  mental  value.  Elsewhere 
there  is  very  little  difference  between  the  resu  Its  calculated  on  thetwo  meshes. 

Fi  gure  19  shows  the  effect  of  the  gri  d  refi  nement  on  the  vel  oci  ty  components  i  n  the  propel  I  er 
pi  ane.  The  resu  Its  are  si  mi  I  ar  to  those  show  n  i  n  Fi  gu  re  12  for  the  effect  of  mesh  refi  nement  on 
the  hexahedral  mesh.  For  the  y  and  z  components  there  is  effectively  no  change  to  the 
simulated  results,  while  for  the  axial  component  the  i ncreased  resolution  has  resulted  in  a 
significant  change  to  the  computed  result  in  the  region  0.0  <  y/  Lpp  <  0.04,  although  this 
change  has  not  significantly  increased  agreement  with  the  experi  mental  result. 


Figure  17.  Effect  of  grid  refinement:  wave  profile  along  the  hull  for  original  tetrahedral  mesh  and 
refined  tetrahedral  mesh 
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x/Lpp 

Figure  18.  Effect  of  grid  refinement:  wave  profile  along  the  line  at  y/Lpp  =  0.172  for  original 
tetrahedral  mesh  and  refined  tetrahedral  mesh 


Figure  19.  Effect  of  grid  refinement:  velocity  components  on  propeller  plane  at  z/Lpp  =  -0.02  for 
original  tetrahedral  mesh  and  refined  tetrahedral  mesh 
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7.  Discussion 

After  having  completed  the  simulation  on  thehexahedral  mesh  it  was  noticed  that  the  hull 
had  been  set  to  an  incorrectvalueoftrimfortheprescribed  speed  of  4.01  knots.  Having  set  the 
hull  to  the  correct  sinkage  at  the  midship  point  it  was  then  rotated  through  an  angle  of  1°, 
rather  than  0.1°,  which  meant  that  the  sinkage  at  the  Aft  Point  was  -0.0096L  and  at  the 
Forw ard  Poi  nt  i  t  w as  +0.0064L .  Th i  s  i  s  si  mi  I  ar  to  the  tri  m  f or  the  sh i  p  travel  I  i  ng  at  a  speed  of 
6.0  knots  rather  than  4.01  knots,  although  the  FP  i  s  hi  gher  i  n  thi  s  case.  Weattempted  to  correct 
this  error  in  the  first  instance  by  going  back  to  the  Gambit  journal  file  and  rotating  the  hull, 
and  the  attached  mesh  near  to  the  surface  of  the  hull  (which  had  initially  required  the  most 
work),  but  thi  s  led  to  substantial  changes  i  n  the  edge  and  face  numberi  ng  scheme  so  that  the 
journal  file  failed  to  run  correctly.  We  believe  that  this  problem  stems  from  a  degree  of 
i  ncompati  bility  bdween  the  Gambit  software  and  the  way  theface  geometry  was  descri  bed  i  n 
the  initial  geometry  file,  which  was  obtained  from  the  Carderock  Division,  Naval  Surface 
Warfare  Center  [29], 

Since  considerable  work  would  have  been  required  to  remesh  the  wireframe  model  from 
scratch  at  the  correct  values  of  sinkage  and  trim  we  decided  to  overcome  this  problem  by 
rotating  the  entire  mesh  relati  veto  the  coordinate  system  (and  therefore  the  inflow  boundary 
conditi  on)  so  that  the  model  was  pi  aced  at  the  correct  si  nkage  and  tri  m.  Thi  s  meant  that  the 
model  was  now  at  the  correct  attitude  with  respect  to  the  flow,  but  this  approach  created 
several  new  problems.  Because  of  the  large  extent  of  the  domain,  and  the  relatively  small 
depth  of  the  ai  r  I  ayer  over  the  sea  surface,  when  the  enti  re  mesh  was  rotated  it  meant  that  the 
outflow  boundary  was  filled  entirely  with  only  onephase(water),  rather  than  the  correct  ratio 
of  thetwo  phases.  This  situation  was  not  allowed  by  theOpen  Channel  boundary  condition. 
The  probl  em  was  easi  I  y  sol  ved  by  si  gni  fi  cantl  y  i  ncreasi  ng  the  hei  ght  of  the  ai  r  I  ayer  u  nti  I  the 
correct  ratio  of  thetwo  phases  was  again  approached.  Another  problem  which  emerged  from 
this  approach  w  as  that  the  hexahed  ral  cel  I  s  al  ong  the  sea  su  rface  were  now  no  I  onger  al  i  gned 
with  theflow.  Although  the  degree  of  misalignment  was  small,  the  effect  on  someofthehigh 
aspect  ratio  cells  was  sufficient  to  cause  considerable  problems  with  the  smooth  running  of 
thesi  mu  I  ati  on.  The  relaxation  factors  had  to  beconsiderably  reduced  for  example  in  order  to 
avoid  frequent  code  crashes. 

The  main  problem  encountered  with  this  approach  however  was  that  the  air/ sea  interface 
became  consi  d  erabl  y  d  i  sru  pted  d  u  ri  ng  the  cou  rse  of  the  si  mu  I  ati  on .  L  arge  amou  nts  of  w  ater 
were  found  to  be  suspended  in  the  air  in  completely  irregular  patterns.  A  remedy  for  this 
problem  was  eventually  found  after  discussions  with  a  Fluent  engineer  [30]  and  consisted  of 
"repatching"  the  appropriate  part  of  the  domain  with  theair  phase  during  the  course  of  the 
simulation.  This  "ad-hoc"  fix  worked  well  on  a  Fluent  machine,  but  was  not  100%  successful 
for  our  simulation,  which  was  run  on  an  eight  processor  Beowulf  cluster  using  Opteron 
processors.  In  particular,  the  problem  was  found  to  persist  on  the  symmetry  plane  no  matter 
how  many  ti  mes  we "  repatched "  the  ai  r  domai  n .  H  owever,  w  hen  the  si  mu  I  ati  on  was  moved 
to  a  quad  core  Intel  processor  this  problem  did  not  occur,  although  problems  still  persisted  on 
the  symmetry  plane. 
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Although  wecontinued  this  simulation  for  a  considerable  length  of  time  thesea  surface  never 
achieved  the  correct  wave  pattern.  Figure  20  shows  that  the  disturbances  on  the  symmetry 
plane  lead  to  the  formation  of  a  non-physical  wave  some  distance  ahead  of  the  bow  of  the 
ship.  This  removed  energy  from  the  real  leading  bow  wave,  which  was  then  considerably 
reduced  in  amplitude.  The  reason  for  this  unphysical  wavesomedi  stance  in  front  of  thebow 
can  be  seen  in  Figure  21,  which  showstheair/  sea  interface  in  the  region  forward  of  the  bow, 
on  which  is  superimposed  I  ines  representing  thepartiti  oning  of  thedomain  amongst  someof 
theeight  processors.  The  occurrence  of  the  unphysical  waveatthe  interface  between  the  two 
parti  ti  on  regi  ons  I  eads  us  to  bel  i  evethat  there  may  be  some  probl  ems  w  i  th  the  paral  I  el  versi  on 
of  the  Fluent  software.  Itshould  however  be  pointed  out  that  this  problem  did  not  occur  when 
the  same  simulation  was  run  by  a  Fluent  engineer  on  a  different  computer  system.  The 
reasons  for  these  discrepancies  are  still  under  investigation  by  Fluent. 

G  i  ven  that  the  si  mu  I  ati  on  resu  I  ts  p  resented  i  n  Secti  on  5  w  ere  ru  n  at  an  i  ncorrect  val  ue  of  tri  m 
it  is  important  to  consider  any  effect  this  may  have  on  the  validity  of  the  results.  Yang  and 
Lohner  [31]  haveused  an  unstructured  grid-based,  parallel  free-surfaceflow  solver  which  has 
been  extended  to  simulatesinkageand  trim  effects  for  the  calculation  of  steady  ship  waves. 
They  applied  their  codeto  the  simulation  of  wave  profiles,  sinkageand  trim,  and  wavedrag 
for  two  hull  forms.  Their  results  show  that  while  achieving  the  correct  sinkageand  trim  is 
i  important  for  accurate  predicti  ons  of  wave  drag,  the  effect  on  the  computed  wave  profi  I  es  i  s 
miniscule.  Given  their  results,  webelievethatthewave  profiles  calculated  inSection  5  would 
not  be  significantly  different  if  thesimul  ati  on  were  to  be  repeated  using  the  correct  val  uefor 
the  trim. 


Figure  20.  Air/sea  interface  showing  the  unphysical  wave  at  some  distance  forward  of  the  bow 
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Figure  21.  Air/sea  interface  showing  the  unphysical  wave  at  some  distance  forward  of  the  bow  with 


superimposed  lines  representing  the  partitioning  of  the  domain  amongst  some  of  the  eight 
processors 


8.  Conclusion 


We  have  used  Fluent  to  calculatethefree-surfaceflow  around  the  DTM  B  5415  naval  ship 
model.  Simulations  were  performed  using  both  a  structured  hexahedral  mesh  and  an 
unstructured  tetrahedral  mesh  of  lower  resolution.  Our  results  show  that  Fluent  is  ableto 
accurately  si  mu  I  ate  the  total  ship  resistance,  near-field  waveshapes,  and  the  velocity  field  in 
the  propeller  plane.  The  total  resistance  coefficient  calculated  on  the  hexahedral  mesh  was 
found  to  agree  with  the  experimental  value  to  within  3.8%.  The  agreement  using  the 
tetrahedral  mesh  was  considerably  less  accurate,  however  this  was  to  be  expected  duetothe 
low  resolution  of  this  mesh  in  the  boundary  layer  region.  The  simulated  waveshapes  along 
thesurfaceofthehull  and  in  the  near  field  computed  on  both  the  hexahedral  and  tetrahedral 
meshes  were  found  to  be  in  good  qualitative  agreement  between  the  simulated  and 
experimental  profiles,  with  thedegreeof  error  from  theresults  on  thetetrahedral  mesh  being 
si  mi  I  ar  to  thosef  rom  the  hexahed  ral  mesh.  Si  mul  ati  ons  of  the  axi  al  vel  oci  ty  component  i  n  the 
propeller  pi  aneshowed  errors  of  approximately  10%ontheoriginal  hexahedral  mesh  and  this 
was  not  significantly  improved  when  the  simulation  was  run  on  the  adapted  mesh.  Good 
agreement  between  the  calculated  and  experimental  results  was  found  for  the  other  velocity 
components.  Similar  results  were  found  for  simulations  on  thetetrahedral  mesh. 

The  above  results  are  encouraging  and  indicate  that  the  Fluent  code  is  a  viabletool  which 
could  be  considered  for  use  in  more  demanding  naval  problems,  such  as  the  computation  of 
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ship  wakes  undergoing  specified  manoeuvres,  for  example.  Use  of  thecodeisfar  from  routine 
how  ever  and  a  n  u  mber  of  si  gn  i  f  i  cant  p  robl  ems  w  ere  encou  ntered  d  u  ri  ng  the  cou  rse  of  th  i  s 
work.  They  are  summarized  here  to  faci  I  itate  future  applications  of  this  codeto  similar  naval 
problems: 

•  TheOpen  Channel  boundary  condition  provides  an  efficient  approach  to  the 
performance  of  these  simulations,  but  it  appears  to  be  incompatible  with  either  a  curved 
domain  boundary  surface  or  with  an  irregular  surface  mesh  over  the  inlet  boundary. 
Simulations  using  theOpen  Channel  boundary  condition  with  a  regularly  mapped 
rectangular  mesh  attheinlet,  which  transitioned  downstream  to  a  less  ideal  mesh 
geometry  via  a  non-conformal  interface,  were  successful. 

•  We  experienced  problems  calculating  an  accurate  sea  surface  shape  when  using  a 
Beowulf  cluster  based  on  theOpteron  processors.  This  problem  did  not  occur  on  aSMP 
machine  using  Xeon  processors.  The  source  of  this  problem  was  not  identified. 

•  Non-physical  behaviour  was  observed  on  the  symmetry  plane  using  both  Opteron  and 
Intel  processors.  These  problems  can  obviously  always  be  avoided  by  removing  the 
symmetry  planefrom  the  simulation  and  integrating  over  thefull  domain.  This  will 
resultin  an  increasein  computational  time  and  resources  however,  so  further  effort 
should  be  madeto  discover  the  source  of  this  non-physical  behaviour. 

•  Fluent  appearsto  be  particularly  sensitive  to  the  angle  of  the  grid  with  respecttothe 
free  surface  water  line  in  simulations  of  this  type.  This  is  far  from  ideal,  but  the  problem 
can  be  avoided  by  remeshing  the  geometry  for  each  individual  sinkageand  trim 
calculation. 
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